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ABSTRACT 

Massive objects sinking within galaxies or dark matter haloes via dynamical friction 
will exchange momentum with central particles, ejecting them from the cusp and re- 
ducing the density of the inner region. We explore parameter space using numerical 
simulations and give empirical relations for the size of the resulting core within struc- 
tures that have different initial cusp slopes. We show that simple energetic arguments 
can be used to predict these scaling laws. As an application we consider the dwarf 
spheroidal galaxy VCC 128 which has a double nucleus separated by less than a hun- 
dred parsecs. If this galaxy has a surrounding cold dark matter halo with central 
density proportional to r~ Y then these objects should sink to the centre of the cusp 
and coalesce in a few million years. We show that the sinking nuclei naturally cre- 
ate a core equal to their current separation if the initial dark matter cusp is slightly 
shallower than a log slope of -0.75 at ~ 0.1% of the virial radius. The sinking objects 
naturally stall at this radius for many Gyrs. This may be indirect observational ev- 
idence for central dark matter cusps shallower than r _1 at the very centres of dark 
matter haloes. 

Key words: cosmology: theory - dark matter - galaxies: dwarf - galaxies: individual: 
VCC 128 - methods: numerical 



1 INTRODUCTION 

Massive objects orbiting within a cuspy mass distri- 
bution are expected to lose mom e ntum and sink via 
dynamical friction ([Chandrasekhai _ 19431: White! 19831: 
Hernquist fc Weinberg 1 19891 : ICapuzzo-Dolcetta fc Vicarl 



20051 ). While spiralling inwards the massive perturber trans- 
fers momentum to central particles/stars etc, moving them 
to a larger orbital radius. One proposed effect of this pro- 
cess is to make a primordially cuspy dark matter dis- 
tribution shallower (|E1-Zant. Shlosman fc Hoffmanl 1200 ll : 
iTonini. Lapi fc Saraccill2006| ). 

In the prevailing ACDM cosmology primordial dark 
matter haloes are cuspy, having an inner density slope 
p(r) oc r 1 with 7 > 1 beyond ~ 1% of the virial ra- 
dius i|Dubinski fc Carlberdll99ll ; lDiemand. Moore fc Stadel 
2005). By contrast, observations of dwarf galaxies seem 
to indicate that they have a cored mass distribu- 



tion l|Sanchez-Salcedo. Reves-Iturbide fc Hernandez! 120061: 
iGoerdt et alj|2006l : iKlevna et al.ll2003lh while controversial 
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evidence for cored mass distributions in dwarf spir al galaxies 
has been debated for over a decade llMooreHl994h . 

A sufficiently compact sinking perturber will transform 
an initially cuspy host halo into a cored one on some scale. 
Once a co re has formed then dynamical friction is no longer 
effective (jGoerdt et al.l l2006h Dynamical arguments show 
that sinking perturber s must stall at the outer edge of a 
core (|Read et al.ll200rj) . 

Here we quantify this stalling behaviour as a function 
of perturber mass M per t and central cusp slope 7. We con- 
sider a much larger range i n M per t and 7 th an in previous 
papers (IGoerdt et al.ll200rj ; [Read et alj|2006t ) and find that 
stalling persists even at very high perturber mass, as found 
also recently bv lGualandris fc Merrittl (|2008l ). We fit empir- 
ical scaling relations using high resolution V-body simula- 
tions and attempt to explain these using simple energetic 
arguments. As an example we attempt to explain the close 
separation of bin ary nuclei within the Virgo cluster dwarf 
galaxy VCC 128 l|Debattista et alj|200rj) . We show that the 
two nuclei are expected to sink into the centre of a cuspy 
halo, and coalesce there on timescales which are extremely 
short compared to the age of the galaxy. We will discuss 
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13 


1.50 
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10.3 
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1.00 


0.009109 


5.01 


15.0 


2.50 x 10 10 
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0.75 


0.017732 


3.91 


19.0 


2.43 x 10 10 


E 


0.50 


0.027746 


3.38 


22.0 


2.40 x 10 10 



Table 1. The parameter list of the five different dark matter 
haloes, we use in our analytical and numerical calculations. 

further applications to multiple nuclei seen in galaxies like 
Andromeda and in cluster cD galaxies in future work. 

This paper is organised as follows: In §2 we describe 
our analytical framework, which is supported using A^-body 
simulations. In §3 we apply our findings to VCC 128 and in 
§4 we give our conclusions. 



2 TRANSFORMING CUSPS TO CORES 

We use the split-po wer law a, /3, 7 model for our initial back - 
ground distribution (lHernauistll990l ; ISahall992l ; IZhacll99r3 ): 



p(r) 



Po 



(r/rsT [l + (r/r s ) a f-< )/o 



(1) 



where po and r s are the central density and scale length 
respectively, 7 is the inner log slope, a is the log slope at 
the 'knee' of the profile, and (5 is the outer log slope. Since 
we focus on the very inner regions of the halo in our analy- 
sis, a and (3 are not critical. We fix them at the commonly 
accepted values for cold dark matter haloes of 1 and 3, re- 
spectively. 

We will take our fiducial model to represent a low mass 
dark matter halo typical of those surrounding dwarf galax- 
ies with a maximum circular velocity u pca k = 50 km s _1 . 
As our application later, we will consider the Virgo clus- 
ter dwarf galaxy VCC 128 which h as an absolute bolomet - 
ric luminosity of Mb = —15.5 mag (|Debattista et al1 i2006). 
We estimate its maximum circular velocity to be w pca k = 
35 - 65 km s -1 usin g the Faber- Jackson relation of dEs 
l|de Riicke et al.ll2005l ). Assuming a concentration of 15 for 
7 = 1.0, which is a com mon value fo r cosmologically moti- 
vated dwarf spheroidals l|Lokasl 120021 ) . we get a r(?; pcak ) of 
10.75 kpc. For our other models we keep ii pca k as well as 
r[v pea k) constant, just varying 7. This leads to the parame- 
ters given in Table [T] and to the circular velocity curves and 
radial density profiles, which are shown in Fig. [1] 

With the above assumptions, VCC 128 is likely dark 
matter dominated at all radii. We estimate its stellar mass 
distribution in two different w ays. Firstly, we u s e the Sersfc 
profile parameters derived in iDebattista et al] j2006) nor- 
malised to give a total luminosity in the B-band of Mb = 
— 15.5 mag: 



E(i?) = Ts/e exp < -b 



GO 



- 1 



(2) 



with b n = 1.9992?! - 0.3271 (jCraham fc Driverl 12005), n = 
0.55, R e = 14.5 arcsec and assumi ng a B-band mass t o ligh t 
ratio Tb = 3 for dE galaxies as in iRead fc Trenthaml (|2005h . 



We then de-project the stellar mass distribution using the 
usual Abel integral equation (that assumes spherical sym- 
metry) : 



p(r) = 

7T 



cE{R) dR 



dR y/R2 _ r 2 ■ 



(3) 



The resulting density and cumulative mass distributions 
are given by the solid lines in Fig. [T] Secondly, we use the 
stellar masses derived from fits to the spectral energy distri- 
bution ( SEP) of the galaxy nucl eus and the whole galaxy as 
given m IDebattista et al.l (j2006l ). These are overlaid on the 
right panel of Fig. [T] (crosses) and give an excellent match to 
the cumulative mass distribution derived from the Sersfc fit 
to the light profile. We assume from here on that the dark 
matter is dynamically dominant and that the stars are to a 
good approximation a massless tracer population. 



2.1 Analytical model 

Previously published work has demonstrated numeri cally 
(|E1-Zant. Shlosman fc HoffmanlliOQll: IRead et"ai1l2006h and 
semi-analytically (jTonini. Lapi fc Saluccill2006t ) that a sink- 
ing massive compact object - a perturber - will transfer en- 
ergy and angular momentum to the background via dynam- 
ical friction, creating a central constant density core from an 
initially cuspy density distribution. Here, we calculate ana- 
lytically how the size of this core will scale for a given initial 
background density distribution and perturber mass. 

The dynamical friction sinking timescale of the per- 
turber can be de rived using the Chan drasekhar dynamical 
friction formula (|Chandrasekharl Il943t ) . Assuming that the 
perturber is always on a circular orbit, and that the back- 
ground particles have a Maxwellian velocity distribution, we 
can calculate the infall rate: 



dr 
dt 



47rlnA(r)p(r)G 2 M p 

crt ' 

v%(r)d [rv c (r)] /dr 



erf 



v c (r) 



2v c {r) 



2-Ka{r) 



exp 



('•) 



2a 2 (r) 



(4) 



where v c (r) is the circular speed at radius r; M pcrt is the 
mass of the in-falling body; In A(r) is the Coulomb logarithm 
[A = bmax/femin] ; p(r) is the density of the dark matter halo 
at radius r according to equation fT}; and a(r) is the one- 
dimensional velocity dispersion of the halo. 

Equation @ assumes that dynamical friction is a lo- 
cal process that proceeds by local momentum exchange 
between the perturber and the background. This cannot 
be the whole story. Firstly, it implies that if p(r) — > 
there will be no dynamical friction. Secondly, it im- 
plies that if p(r) — » const, nothing special should hap- 
pen and dynamical friction should proceed as for any 

other p(r). Both of these statements conflict with numcr- 
1 1 I 1 I 1 

ical results (Lin & Trcmainc 1983; Trcmainc & Weinberg 
19841 : IColpi. Mayer fc Govcrnato 1999; Goerdt et al. [2006 



Read et all 2006). The implication is that dynamical fric- 
tion is not a local process, but rather i s driven by globa l 
resonances (|Tremaine fc Weinberg] Il984l : IRead et al.ll2006l ). 
In the former case, resonances (and therefore friction) are 
present even outside of a galaxy where p(r) = 0. In the 
latter case, p = const, is a harmonic potential which is es- 
pecially resonant. This allows the perturber and the back- 
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Figure 1. Circular velocities (left) radial density profiles (middle) and cumulative mass (right) for haloes A, B, C, D, E. All five of them 
peak at 50 km s — 1 at 10.75 kpc. The horizontal lines on the right panel show the perturber masses M per t we use in the simulations 
( EI2.2I I, The solid lines and the crosses show estimates of the baryonic (stellar) mass distribution of the dSph galaxy VCC 128. 



ground reach a stable state with no net momentum exchange 
( |Kalnaid[l973 : iRead et al.ll200rj| ). 

Despite the above difficulties, in most situations, equa- 
tion @ works remarkably well. This is because the non-local 
effects enter only logarithmically through In A. We use this 
fact in our calculation, below. We assume that the perturber 
transfers all of its energy and angular momentum locally in a 
spherical shell. Provided that we are not in the regime where 
the local density tends to zero, or the background and the 
perturber are especially resonant this will be a good approx- 
imation. 

Equation Q assumes the perturber starts and remains 
on a circular orbit. As a result, the energy transferred to the 
background particles in moving from radius r + 8r to r does 
not depend on the dynamical friction process itself. All we 
need to know is that the perturber moves inwards from one 
circular orbit to the next. In this case, we have: 



SEt, 



_d_ 

dr 



1 9 

-Mpertfc + Mpert$(r) 



5r 



G dM(r) GM(r) + d$(r) 



2r dr 



2r 2 



dr 



Sr (5) 



where M(r) is the mass of background particles enclosed at 
radius r. We assume that this energy is injected into back- 
ground particles contained within a spherical shell about the 
perturber. These gain specific energy AE given by: 



Aitpr SrAE = SEtmns 
and rearranging gives: 
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4irpr 2 



G AM {r) GM(r) d$(r) 



2r dr 



2r 2 
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dr 



(6) 



(7) 



We may solve equation ([7} exactly if we specify the distribu- 
tion of background particles. A particularly useful choice is 
the split power law models given by equation |TJ with a — 1 
and /3 = 4: 



p(r) = 



M t (3- 7 ) 



M(r) = M t 



(0 



r s 

7-3 



7—4 



(8) 
(9) 



$(r) = 



GM t 
r.(2-7) 



(-7) 



7-2 



- 1 



(10) 



Where M t is the total mass, r s is the scale length, and 7 is 
the central log density slope. This gives: 

GM P 



AE 



4-7 + 



(11) 



2(3 -7)r 

We assume that this energy is injected over one dynami- 
cal time and that the system then rearranges itself on a 
timescale longer than this. In this impulse approximation, 
prior to rearrangement, the potential 3>(r) remains constant, 
while the specific kinetic energy of the background increases 
from Ti —* Tf. We can obtain an estimate of 71 if we assume 
that all of the background particles move on circular orbits: 



Ti (r) 



2 



GM t 



2r 

Now, if we write: 
T f (r)-Ti(r) _ 



(0 



1 + 



(12) 



(13) 



Then the radius at which k ~ 1 is the radius at which the 
background specific kinetic energy doubles. 

After the above energy injection, the system must then 
return to virial equilibrium. Assuming no mass escapes dur- 
ing this phase, the specific kinetic energy must be converted 
to potential energy and the system will expand. The maxi- 
mal effect leads to a constant density core since density can- 
not decrease inwards for a stable system. We assume that the 
radius of this core r c is related to the radius at which the 
specific kinetic energy doubles. Combining equations (|11[) , 
()12[) and (|13|) and assuming r c < r s , we obtain for r c : 



Mr, 



■7 1 



Mt 3 - 7 k 



(14) 



Where k, remains as a fitting parameter. We find k = [13.32, 
2.192, 0.441, 0.224, 0.167] provides a good fit for haloes A, 
B, C, D, E to the simulations presented in £12.21 

Note that we have not considered the angular momen- 
tum imparted to the background by the perturber. This may 
be calculated similarly to the above and is important for de- 
termining the final orbit structure l|Tonini. Lapi fc Saluccil 
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2006), but not the above scaling relations. Note also that 
very similar scaling to that in equation (|14[1 may be obtained 
by simply comparing the enclosed mass within r c with the 
mass of the perturber. Taking M pcr t ~ M(r c ) and assuming 
r c <C r s , we hav e r c /r s oc (Mpert/ Mt) 1 ^ 3-7 - 1 . However, as 
demonstrated in iRead et all (|2006D . the stalling behaviour 
that occurs when p ~ const, is not simply the result of 
M peT t ~ M(r c ). For large cores, stalling occurs even when 
this criteria is not met. Finally, note that r c does not depend 
on the initial radius of the perturber, only on M pC rt and the 
distribution of the background. 

Our result agrees well with previous studies. 



lEl-Zant. Shlosman fc Hoffman! (|200lf ) use a semiana- 
lytical Monte Carlo approach based on the Chandrasekhar 
approximation to estimate core sizes. They let 100 — 500 
perturbers, which have a combined mass of 10% of the host 
halo, sink into a 7 = 1.0 halo. They find r c « r s . Putting 
Mpcrt /Mt =0.1 into equation (fT4)l we find, using k = .441 



for a 7 = 1.0 halo, r c = 0.58r s . iGoerdt et alT(|2008l ') let 
10 live globular clusters, which have a combined mass of 
0.28% of the host halo, sink into a 7 = 1.0 halo. They do 
not measure the core size explicitly but in their Fig. 1 the 
initial dark and the final dark lines depart at around 200 
pc. Putting Mp er t/A/t = 0.0028 into equation (JT4} we find, 
using k — 0.441, r c = 150 pc. Given the uncertainty on «, 
our estimates agree very well with what has already been 
published in the literature. 



2.2 N-body simulations 

In this section, we perform iV-body simulations with the 
haloes A, B, C, D, E presented in Table [1] to study 
the effects of different mass perturbers on modifying the 
density profiles. Equilibrium TV-body representations of 
th ese haloes were created with th e algo rithm described 
in iKazantzidis. Maeorrian fc Mqorel (I2004TI using the mul- 



timass technique of IZemp et al. l|2008h .~This gives struc- 



tures that have an effective resolution in the region of in- 
terest equivalent to using ~ 10 10 single mass particles. We 
ran a grid of simulations with M pcrt = [10 5 , 5 x 10 5 , 2.5 x 
10 6 , 10 7 , 5x 10 7 ] M for each of the haloes A, B, C, D, E cor- 
responding to varying the initial central log density slope: 
7 = [1.75, 1.5, 1, 0.75, 0.5] (see Table [TJ. The perturbers 
were typically started at a radius of 0.4 kpc, except for the 
heaviest perturbers which we had to start further out for 
reasons which will be explained in ^2.2.11 (We could have 
started all nuclei at the outermost radius and obtained the 
same results, but this would have produced avoidable com- 
putational costs. See Fig. [5]) All simulations are shown us- 
ing circular orbits but similar scaling laws were found using 
more eccentric orbits. The resulting trajectories are given in 
Figure [2] and the resulting stalling radii in Figure [3] 



2.2.1 Kickback and stall 

For all trajectories there is a kickback which occurs after a 
first point of closest approach (fpea). The perturber moves 
away for a while, reaches a maximum, and then returns to 
a second point of closest approach (spca), where it finally 
stalls. The distance of the perturber from the centre of the 
halo is the same up to the noise level at fpea and at spca. We 




Figure 5. Simulated position of a single perturber with identical 
mass (5 X 10 7 Mq) but different initial radii within halo D as a 
function of time. The stalling radius is indicated by the solid line. 



define the height of the kickback to be the difference between 
the lowest and the highest point during the kickback period. 
It is interesting that the local minima to the right and to 
the left of the kickback, fpea and spca, are approximately 
equal in height. 

The kickback occurs at ~ r c because M pcrt ~ M(r c ) 
(see £|2.1[1 . At this point, the acceleration on the perturber 
due to the background is equal to the acceleration on the 
background due to the perturber, and the centre of mass of 
the system is significantly displaced. The background, unlike 
the perturber, however is not solid and rapidly rearranges 
itself to form the central constant density core, after which 
the perturber stalls. We assume that the duration of the 
kickback is proportional to the crossing time within the core 
^kickback = X*cro SS , which gives: 



tkickback = X~ 



(15) 



*«c (r c ) ' 

Using equations (|9}, (|14[) and assuming r c < r B , this gives 



tkickback = X 



GM t 



M p 4- 



7 1 



M t 3 - 7 k 



The height of the kickback may then be estimated as: 

/ikickback = tkickback ^c (?"c) = 4>T C 



(16) 



(17) 



The duration and height of the kickback as a function of 
Mpcrt are given in Fig.rj. We find X = 80.56 and = 0.5587 
provide good fits. 

In Fig. [S] we show the stalling radius r c as a function 
of the initial radius n. As expected from the analytic calcu- 
lation in ij2.ll r c does not depend on n as long as n > r c . 
Only the simulation which starts with n < r c stalls well in- 
side r c . Since all of the simulations had to be started outside 
of r c , the very heavy perturbers had to be started further 
out than the light perturbers. This does not affect any of the 
results, presented here. Both the height and the duration of 
the kickback are independent of the initial radius as well, in 
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Figure 3. The stalling radii as a function of M per t for varying 7 (left) and as a function of 7 for varying Af pcr t (right). Analytic models 
are also shown. 



the sense that there is no monotonic relation between n and 
either the height or the duration of the kickback. 



2.2.2 Core creation 

The density profiles of the respective host halo at fpea and 
at spca are plotted in Fig. [6] and Fig. [7] One can clearly 
see that the density distribution changes significantly. The 
haloes settle down into a new equilibrium between fpea and 
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Figure 4. Duration (left) and height (right) of the kickback as a function of M per t for haloes B, C, D, E, together with our analytic 
models. We do not show halo A on these plots because the data were too noisy to give reliable data points. 



spca. The radial density profiles change from being cuspy 
to having a core: as the perturber mass increases the initial 
cusp is gradually flattened to a systematically higher de- 
gree until a constant density core is attained. Note that our 
previous analytic estimates give an estimate as to the radius 
within which the density profile is affected by the perturber. 

The change from cuspy to cored central density is im- 
portant for the expected annihilation signal from a weakly 
interactive massive (WIMP) dark matter particle, since 
the annihilation signal goes as the density squared. The 
net flux coming from WIMP annihila tion is given by (e.g, 
lKoushiappa3l2006l : iGoerdt et all I2007T I : 



F 



4nr p(r) dr 



(18) 



where the dependence of the flux on the WIMP mass and 
interaction cross section is wrapped up inside the constant 
k. The lower bound r m i n is defined as the central region of 
the host halo in which the neutralinos have already annihi- 
lated. The required number density for this to happen can 
be estimated using: 



th 



1 

nav 



(19) 



where th « 13 Gyrs is the Hubble time, av i 
is a typical cross section and n is the number density of 
neutra linos. For more details see lCalcaneo-Roldan fc Moord 
(1999). The minimum radius can now be computed (for a 
sinking perturber of a given mass) by comparing this mini- 
mum number density with the density profile in Fig. [7] As- 
suming a WIMP mass of 100 GeV and deploying the above 
mentioned density profile, r m j n is of order 10 -14 pc. Fig. [8] 
shows the resulting annihilation flux. It is more or less in- 
dependent of the assumed WIMP mass, because this mass 
only goes into the calculation of r m j n , which is very small. 
For typical core sizes of r c ~ 0.1r s , core creation can lead 
to a decrease in flux of up to a third. A much weaker effect 
would result from a single sinking star; for example a WMq 



0.9 



0.8 



0.7 



0.6 



0.5 



0.4 



0.3 



Y= 1.25 
y= 1.00 
y=0.75 
y=0.50 



0.05 



0.1 



0.15 



0.2 



, /r. 



Figure 8. The flux of annihilation products from different host 
haloes after the central cusp core transformation relative to the 
untransformed initial cusp. Typical core sizes (r c fa 0.1r s ) can 
lead to a decrease in flux of up to one third. The solid line shows 
7 = 1.25 and is not related to any of the haloes used through- 
out the rest of the paper (the corresponding values for 7 ^ 1.5 
diverge) . 



star would creating a core of radius 0.34pc in our fiducial 
halo. 



3 APPLICATION TO VCC 128 

In this section, we apply our results to VCC 128, a dwarf 
spheroidal galaxy (dSph) at the outskirts of the Virgo clus- 
ter. A very close binary nucle us has recently been discovered 
there (jPebattista et al1l2006l ). The two nuclei are very simi- 
lar in their appearance with masses estimated to be around 
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Figure 6. Density profiles of the host halo at fpca for the different nuclei masses M port . From left to right upper panels show haloes A, 
B, C and lower panels haloes D and E. 
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Figure 7. Density profiles of the host halo at spca for the different nuclei masses M per t- From left to right upper panels show haloes A, 
B, C and lower panels haloes D and E. 
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5 x 10 5 Mq . These authors found altogether 3 dSph galaxies 
with double nuclei, looking at HS T archival images of 50 
dwarf spheroidals from a survey bv lLotz. Miller fc Ferguson! 
|2004 ). T he projected distance of the two nuclei in VCC 128 
is 32 pc. iDebattista et all |2006Tl suggest, because the two 
nuclei have very similar colours and magnitudes that this 
could be evidence for a nuclear disc around a supermas- 
sive black hole ( SMBH), a situation as in N GC4486B 
dLauer et al.lll996T ) and similar to the one in M31 (|Tremainel 
1 19951 ). However, it is not confirmed t hat such an SMB H can 
exist in a dwarf galaxy like VCC 128. lFerraresej (|2002h found 
from a sample of 36 galaxies, tentative evidence that SMBH 
formation becomes inefficient in haloes below a dynamical 
mass of ~ 5 x 10 n Ma, though more recent work may 
sugge st otherwise l|Ferrarese et al.l 120061 ; IWehner fc Harrisl 
2006). As such, it is interesting to explore alternative expla- 
nations. 

Figure shows the time taken for two nuclei that have 
masses 5 x 10 5 M© and initial separation of 44 pc to sink 
via dynamical friction in halos B,C,D and E, assuming that 
equation U fully describes the friction process. In all haloes 
with cusp slopes shallower than 0.75, the nuclei coalesce 
rapidly. This suggests that, were equation [4] the whole story, 
it would be very unlikely to observe 3 close nuclei out of 50 
dSph galaxies in Virgo. However, as we have demonstrated 
in the previous section, the merging nuclei will create a small 
core and stall indefinitely, leading to a much higher proba- 
bility of observing double nuclei. We quantify this in detail 
below. 

We assume that all dwarf galaxies for which no double 
nuclei have been observed: (a) do not have double nuclei at 
their centres; (b) lost their double nuclei because of a sink-in 
process due to dynamical friction, and that all dwarf galax- 
ies, for which double nuclei actually have been observed; 
and (c) have double nuclei as close as the ones in VCC 128 
(this assumption is made not only for the sake of numerical 
convenience, but also because of lack of precise observational 
measurements) . It is important to note that there are uncer- 
tainties in these assumptions: Out of the 47 dwarf ellipticals, 
which do not have double nuclei, there might be some which 
actually do have such nuclei but the latter are too close to 
each other to be observed. Moreover, not all dwarf galaxies 
necessarily developed double nuclei during their lifetime. 

Throughout this section we will assume that out of two 
orbiting nuclei one sinks to the centre of the underlying dark 
matter distribution and the other one is orbiting it after- 
wards. So we interpret the distance between the two nuclei 
as the distance of one nucleus to the centre of the dark mat- 
ter distribution. We performed several iV-body simulations 
to verify that this assumption holds. 

Using equation Q, we can calculate the time it would 
take for one of the nuclei to sink to the centre of the system. 
In all our calculations we assume a constant Coulomb log- 
arithm In A(r) = 4.0 (see iPenarrubia. Just fc Kroupa|[2004l . 
for a detailed discussion of In A varying with r). We as- 
sume that the vector connecting the two nuclei has a ran- 
dom orientation. We define inclination in a way that it is 
equal to 0° if there is a right angle between the line of sight 
and this vector. The probability density of its inclination is 
^(i) — cos(i). From this, it follows that in P ang i c = 70% of 
the cases the inclination is lower than 43° or the de-projected 
distance between the nuclei is lower than 44 pc. To make a 



halo B 
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halo E 
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Figure 9. Analytical computed position as a function of time, 
during the Chandrasekhar sink-in period for four of our haloes, 
assuming an initial true distance of the nucleus from the centre 
of 44 pc. 



start, we will use this value as the true distance. For halo B 
the total sink-in time is computed to be 7.4 Myr, Halo C's 
sink-in-time is 13.9 Myr, Halo D's sink-in time is 38.1 Myr 
and Halo E's sink-in time is 171 Myr. The sink-in behaviour 
of the nuclei for all four cases can be seen in Fig. [9] We as- 
sume that the current age of the galaxy is 10.0 Gyr. So the 
time it takes for the two nuclei to coalesce is about 0.075%, 
0.14%, 0.38% and 1.7% of the lifetime of the galaxy, respec- 
tively. This means that it would be highly improbable to 
observe binary nuclei in the present. We want to elaborate 
a bit further on the comparison of these probabilities. The 
sample we are working with has 50 galaxies in total, out of 
which 3 have double nuclei. When we look at one random 
galaxy at a random time for halo D, the probability of seeing 
a double nucleus is Pdn = 0.0038. The probability of seeing 
more than 2 galaxies with double nuclei in a sample of 50 is 
given by the cumulative binomial distribution: 



Pn 



d„>2 



£ 

j=o 



(Pd„) J (l-Pdn) 



50-j 



(20) 



For the values of our halo D we have P ndn >2 = 0.0941% 
or alternatively P ndn <3 = 99.91%. This is valid under the 
assumption, we made earlier, that the inclination angle i is 
smaller than 43° . But this assumption holds only in P ang ic = 
70% of all the cases. Therefore we can exclude halo D with 



a probability of P. 



angle 



= P 



angle-* n dn <3 



69.9%. 



In order to obtain exclusion probabilities, which are in- 
dependent of an arbitrarily chosen inclination angle i, one 
has to integrate the i dependent probability of seeing less 
than 3 galaxies with double nuclei in the sample Pn dn <3(i) 
over all angles normalised by their probability density: 

/■tt/2 

Pcxdudc = / cos(i)P ndn<3 (i)di (21) 
Jo 

Table [2] shows those values for the different haloes. A clear 
trend becomes apparent. The excludability increases with 
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Halo 



B 


1.50 


99.31% 


A 


2.70o 


C 


1.00 


98.89% 


A 


2.54ct 


D 


0.75 


97.99% 


A 


2.32ct 


E 


0.50 


90.63% 


A 


1.68o 



Table 2. Statistical implications of neglecting the cusp-core 
transformation mechanism explained in [[2] for the different halo 
models derived by our analytical calculations. 
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Figure 10. Separation of the two perturbers as a function of time, 
assuming an initial true separation of 800 pc and perpendicular 
orbits. The horizontal line shows the observed projected distances 
in VCC 128. 



increasing 7, as we might expect. So there is an upper limit 
on 7. Looking carefully at the probabilities given in Table [2] 
we must come to the conclusion that in absence of the cusp- 
core transforming mechanism described in the last section 
the inner log density slope of VCC 128 must be extremely 
shallow. 

We ran an additional series of simulations with two per- 
turbers on either coplanar or perpendicular orbits. Again we 
used the haloes presented in Table [T] The stalling behaviour 
is shown in Fig. 1101 In haloes B, C and D, the nuclei do not 
stall. Only in halo E do they stall above 32 pc. These results 
suggest for slopes of 1.0 or steeper, this stalling mechanism 
does not crucially change the results we derived analytically 
in this section. However for slopes shallower than 1.0 it af- 
fects these results quite dramatically. 

To summarise our findings: We can indeed exclude all 
steep inner log density slopes of 0.75 or higher. The most 
likely explanation for the observation is however that the 
two nuclei transformed an initially steeper density profile 
with an inner log density slope between 0.75 and 0.5 into a 
core during its initial sinking period. The nuclei would then 
stall at the edge of this freshly created core. A very shal- 
low inner log density slope like 0.5 or smaller seems to be 
quite unlikely as well, because the globular clusters would 



have stalled far further away then 32 pc (at 100 - 200 pc; 
see Fig. I10[) . Therefore one can deduce that an underlying 
dark matter halo with 7 ~ 0.75 — 0.5 at ~ 1% of the virial 
radius provides the best explanation for the shown obser- 
vations. One should note here that this cusp-core transfor- 
mation mechanism does not explain the dynamical friction 
timescale prob lem of Fornax within a ACDM (7 = 1.0) halo 
(see Fig. 3 of iGoerdt etafl £2006)), because in that case 
the globular clusters in Fornax are far too light and far too 
distant from the centre of their host galaxy for this cusp- 
core transformation mechanism to play a significant role. 
Equation (|14jl predicts for the halo they are using, assum- 
ing k = 0.441 for a 7 = 1.0 halo, a core of only 44 pc. This is 
far smaller then the radius of the innermost globular cluster 
in Fornax, which is at 240 pc. 



4 CONCLUSIONS 

A sufficiently compact sinking perturber will transform an 
initially cuspy host halo into a cored one on some scale. 
Once a co re has formed then dynamical friction is no longer 
effective (IGoerdt et al.ll2006l ). Dynamical arguments show 
that sinking perturber s must stall at the outer edge of a 
core l|Read et al.ll2006l ). 

In this paper we quantify this stalling behaviour as 
a function of perturber mass and central cusp slope. We 
give empirical scaling relations using high resolution iV-body 
simulations and attempt to explain these using simple ener- 
getic arguments. We describe a mechanism which is capable 
of transforming primordial cusps into cores. At first sight, 
this mechanism looks promising for solving the cusp core 
problem because it naturally explains cores in a ACDM uni- 
verse. From equation (|14|l . a ~lkpc sized core will form 
from a single perturber having 1% of the mass of its host 
(assuming k = 0.441 for 7 = 1.0). This essentially r e covers 
the earlier results of El-Zant. Shlosman fc Hoffman! (|200lh 
and lEl-Zant et"aL I <|2004h . However, we require a point mass 
with 1% of the host mass - that is, for example, a 10 7 Mq 
point mass in a 10 9 Mq dwarf galaxy - which seems unreal- 
istic. But there is a further problem. The effect is strongest 
only when the mass of the perturber matches the enclosed 
mass of the host galaxy. Thus, unless the perturber is a black 
hole, a bright dense nucleus will form making the dark mat- 
ter dynamically irrelevant. Such galaxies undoubtedly ex- 
ist, but they are not the typical L SB galaxies fo r which the 
cusp-core problem is most famous (de Blok 2005). All of this 
explains why detailed numerical simulations that model the 
dissipation of gas and star formation do not see cusps turn 
into cores, de spite the stars bein g over-concentrated due to 
over-cooling (iMaccio et alj|2006h . 

Although not necessarily important for the cusp-core 
problem, the cusp-core transformation mechanism we de- 
scribe in this paper could explain the abundance (3/50) of 
dwarf spheroidal galaxies (dSphs) in Virgo with close double 
nuclei. We focused on the dSph VCC 128 which has a dou- 
ble n ucleus of projected separation 32 pc (|Debattista et all 
2006). We suggested a new model where the double nucleus 
is formed by two infalling globular clusters that create a 
smaller constant density core and then stall. We demon- 
strated that VCC 128 is likely dark matter dominated and 
placed constraints on the distribution of dark matter re- 
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quired for this core-stalling mechanism to be effective. We 
found that: 

• The inner log density slope 7 of the dark matter halo 
of VCC 128 must be less than 0.75 at ~ 0.1% of the virial 
radius. For 7 > 0.75 initially, the dynamical friction sink-in 
time is so small in comparison to the lifetime of the galaxy 
that we run into a fine tuning problem. For 7 < 0.75 initially, 
the nuclei create a central constant density core of separation 
~ 40 pc after which they stall. 

• The initial inner log density slope 7 must be greater 
than 0.5. Otherwise the globular clusters would stall well 
beyond their current projected separation of 32 pc. 

The above may be indirect observational evidence for 
central dark matter cusp slopes shallower than r _1 near the 
centres of dark matter haloes. Unfortunately, published nu- 
merical simulations of CDM halo structure have not quite 
reached this resolution (|Graham et al.| [2006). which requires 
following over a billion particles within the final halo. 
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